Energy gap tuning in graphene on hexagonal boron nitride bilayer system 
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We use a tight binding approach and density functional theory calculations to study the band 
structure of graphene/hexagonal boron nitride bilayer system in the most stable configuration. We 
show that an electric field applied in the direction perpendicular to the layers significantly modifies 
the electronic structure of the whole system, including shifts, anticrossing and other deformations 
of bands, which can allow to control the value of the energy gap. It is shown that band structure 
of biased system may be tailored for specific requirements of nanoelectronics applications. The 
carriers' mobilities are expected to be higher than in the bilayer graphene devices. 
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I. INTRODUCTION 

Graphene, a one atom thick graphitic carbon, since 
its fabrication in 2004i has attracted a huge attention 
due to its exceptional electronic properties^ and currently 
is a subject of extensive research in condensed matter 
physics. Due to extremely high carrier mobility graphene 
is considered to be an ideal material for carbon-based na- 
noelectronics. The production of large enough graphene 
sheets makes it a promising candidate for direct appli- 
cations, especially as a channel material in field effect 
transistors (FETs). 

Achieving large on/off current ratios is essential for 
logic devices. It has been demonstrated that semimetal- 
lic single layer graphene (SLG) is not a suitable choice for 
FETs constructions^ since only an on/off current ratio of 
about 4 has been observed.^ Although graphene nanorib- 
bons field effect transistors (GNRFETs) have been pro- 
duced, exhibiting excellent properties^i^, mass produc- 
tion of GNRs-based devices is beyond abilities of current 
litography techniques J. 

The opening and external tuning of energy gap be- 
tween valence and conduction bands in Bernal stack- 
ing bilayer graphene (BLG)^i^ reveals that it can have 
the largest potential for logic applications. It has been 
demonstrated experimentally using angle-resolved pho- 
toemission spectroscopy (ARPES) that energy gap tun- 
ing can be achieved in epitaxial graphene.-iS Charge 
transfer from SiC substrate and additional doping by 
potassium atoms on top of the graphene sample in- 
duce effective external field controlled by changes of 
dopants amount. However, recent experimental studies 
on very large number of samples produced in an inte- 
grated technologjsii show that switching off the conduc- 
tion to desirable level is still a very difficult task in epi- 
taxial graphene. 

Energy gap tuning by electric field effect has been 
also detected in exfoliated bilayer graphene deposited on 
Si02/Si wafer controlled by using double gate device i^^d^ 
The optical measurment o^^'^^ have confirmed predictions 
from tight binding calculations that energy gap saturates 



and does not achieve more than 71 ~ 0.35 eV.^^ A re- 
cent study of carrier transport in state-of-the-art bilayer 
graphene FET device^ indicates that a transport band 
gap of >130 meV and on/off current ratio of around 100 
at room temperature is observed. 

In bilayer graphene, in contrast to single layer sheet, 
the Dirac fermions are massive, with an effective mass 
m* ~ 0.03 TTic (where mo is the bare electron mass);^ 
Since the exceptional conical dispersion is not then pre- 
served, an important question arises whether it is possible 
to open a band gap and still keep the extraordinary elec- 
tronic properties of SLG. Several proposals of band gap 
engineering have been recently given, including surface 
bonding (graphane), isoelectronic co-doping with boron 
and nitrogen atoms and alternating electrical or chemi- 
cal environment (by using Li and F atoms) It has been 
demonstrated, however, that these systems may not offer 
higher carriers' mobilities than BLG-based semiconduc- 
tors. We believe that the crucial issue of band gap tuning 
and at the same time preserving the exceptional linear 
dispersion can be achieved by using suitable substrates. 

We reintroduce here the idea of graphene on hexago- 
nal boron nitride (h-BN) given in Ref. [l^. In contrast 
to graphene, the isoelectronic h-BN crystal with strong, 
partially ionic sp^ hybrydized in-plane bonding between 
boron and nitrogen atoms, has a wide energy gap re- 
sulting from a significant difference of electronegativities 
between onsite atoms. It emerges in a natural way that 
structural and band structures similarities between these 
two materials can lead to amazing physical properties of 
a combined system. 

It has been showrti^ that the most stable configura- 
tion of a graphene on h-BN substrate has one carbon 
atom on top of a boron atom while the other one sits 
above a center of boron nitride ring. The density func- 
tional theory (DFT) calculations indicated the opening 
of a gap of 53 meV at the Dirac points resulting from the 
graphene-substrate interaction. On the other hand, as in- 
dicated in Ref. Hi], one can calculate the band structures 
of graphene and hexagonal boron nitride diagonalizing 
a hamiltonian of nearly the same form. This similar- 
ity allows to calculate conduction and valence bands of 



2 



graphene/h-BN bilayer system in almost the same man- 
ner as for graphite^^ BLG and FLG.<22"— The extension 
to the case of boron nitride sheet in place of carbon bot- 
tom layer is straightforward and perpendicular external 
fields can be easily taken into account. 

In the present paper the electronic properties of a bi- 
ased graphene/h-BN system are studied within a tight 
binding (TB) approach verified by ab initio calculations. 
The band structure is calculated both in the presence 
and in the absence of a perpendicular electric field. We 
study how the gap can be externally tuned by the elec- 
tric field effect. The advantage of the present approach 
is that substrate, directly incorporated into the model, is 
an integral part of system which is indispensable in mod- 
elling the band structure. It is shown that graphene on 
h-BN layer can be an alternative system with externally 
tunable energy gap. 

The paper is organized as follows. In Sec. II we dis- 
cuss the band structure evolution from bilayer graphene 
towards the graphene on hexagonal boron nitride in the 
framework of the nearest-neighbour tight binding formal- 
ism. Since the hopping parameters needed in the tight 
binding description are taken from experimental or ab 
initio calculations data (TB-LDA^ TB-GW^^), the fit 
of TB energy bands to DFT results is needed and pre- 
sented in Sec. III. The final remarks and perspectives are 
reported in Sec. IV. 



II. TIGHT BINDING CALCULATIONS 

In periodic table of elements carbon is placed between 
boron and nitrogen, thus they differ only in the num- 
ber of 2p electrons. Since the honeycomb structure of 
graphite originates from sp^ hybridized in-plane bond- 
ing, the proximity of electron energies of 2s and 2p states 
in carbon, boron and nitrogen^ results in formation of 
similar structure in the case of boron nitride. Graphene 
and h-BN have identical lattices, but, in contrast to zero 
energy gap in graphene, the gap in h-BN is of the order 
of 5 eV depending on the number of layers. The differ- 
ence between two materials can be explained within TB 
approach. 

Both crystals have two dimensional hexagonal struc- 
ture shown in Fig. [T](a). A unit cell contains two sites, 
one of type A and the second one of type B. The linear 
combinations of lattice vectors: 
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generate only A sites while sites in B sublattice are gen- 
erated by nifli -|- 712(12 + (5, where d has to be choosen as 
one of the three nearest neighbours vectors: 



The quantity a stands for a lattice constant. The corre- 
sponding Brillouin zone (see Fig. [l](b)) is a hexagon in 
reciprocal space with basis vectors given by: 
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Two inequivalent corners of the Brillouin zone K and K' 
can be chosen as follows: 
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Although there is a slight difference between lattice con- 
stants of graphene and hexagonal boron nitride, we take 
the same value for both materialsfi^ 

We present a simple model of graphene on hexagonal 
boron nitride in Fig. [l](c). The top layer consists of car- 
bon atoms, while the bottom layer is a hexagonal boron 
nitride sheet. As indicated in Ref. [l^, it is the lowest en- 
ergy orientation of the graphene layer with respect to the 
h-BN with the equilibrium distance between the sheets 
c = 3.22 A. A two dimensional Brillouin zone is shown 
in Fig. 1 (b). One can regain the AB stacked bilayer 
graphene replacing boron and nitrogen by carbon atoms 
and substituting interlayer distance by c = 3.35A. 

Due to sp^ hybrydization the resulted three a bonds 
are located in the graphene plane while two overlapping 
2pz orbitals form tt bond perpendicular to the sheet. 
Therefore a and tt orbitals can be analyzed indepen- 
dently. We restrict our considerations only to tt energy 
bands because they are responsible for unique proper- 
ties of graphene. Two Bloch functions, constructed from 
atomic orbitals for two inequivalent carbon atoms A and 
B, provide the basis functions for single layer graphene. 
For bilayer system we assume a basis as 2pz orbitals of 
respective atoms with Bloch factors: 
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where Rp-^^ and Ra2, ^'^^ the positions of in- 

equivalent atoms in the first and second layer, while (pa-^ , 
and (f)ci2 J (t'P2 "^Pz orbitals in respective atoms. 
We use atomic valence orbitals and orbitals with close 
values of energies^ ([5]) to construct the eigenfunctions: 



k k ^ ' 



k k ^ ' 



(6) 



«(- 



2V3' 2^ 



(2) 



3 



32 




FIG. 1. (Color online) (a) Lattice structure of graphene (or h-BN) constructed by two interpenetrating triangular sublattices, 
(b) corresponding Brillouin zone and (c) lattice structure of graphene (black points) on h-BN. Boron atoms are denoted as red 
(darker) circles, nitrogen atoms as blue (lighter) ones. 



According to the procedure given in Ref. [2^ we intro- 
duce the following quantities representing the energies of 
electrons of 2p level in the crystal: 
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and in-plane and an out-of-plane hopping parameters: 
7o = / dr {r - R^,)H<t>p, (r - - 6) (8) 
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and we define the overlapping matrix in the form: 
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where so, Sq and s\ denote overlapping integrals which 
are assumed to be small. In fact, since we study the band 
structures in the vicinity of the Brillouin zone corners, we 
can neglect all of them.— 

The eigenvalue problem is simple in case of monolayer 
system and we start with providing a short outline of the 
obtained results. 



71 = y dr i?aji?0&(r - R^, - c) (10) 

where iJ is a hamiltonian of the system, c = c Cz and 
only the nearest neighbours interactions are assumed in 
the system. We construct then a hamiltonian in the basis 
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Graphene and boron nitride monolayers - an 
outline 



In case of graphene monolayer the hamiltonian (jlip 
reduces to right bottom block with the orbital energies 



of the 2pz level in carbon lattice t^-p^ ^ ^2p, — ^ 
have the solutions: 



i?±(fc) = 6^p^±7o|/(fc)| 
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where ± signs define two branches of the energy disper- 
sion curves and E-{k), E^{k) are called bonding tt and 
antibonding tt* energy bands, respectively. It is seen from 
definition that at the K and K' points it holds: 
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and a zero band gap is easily recognized. It means that 
the electrons in two sublattices do not interfere if they 
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are at the Brillouin zone corners and it is a direct conse- 
quence of the symmetry of the system. It is worthwhile 
to remind that this does not change if hopping between 
non nearest neighbours is included in the model. 

The hamiltonian (fTTj) reduces to left top block for the 
case of boron nitride monolayer. The solutions can be 
written in a simple form: 

E±{k) = I (e^,^ + 6^ ± (16) 
with 

defining the band gap at K point. Then, in the corners 
of the Brillouin zone we have: 



In this case the symmetry between A and B sublattices 
is broken because the energies of the onsite atoms are not 
the same representing boron and nitrogen, respectively. 
The band structure is similar, but at the K {K') point 
a gap with the magnitude Eg opens. In the first approx- 
imation, it emerges that energies e2p, and ejp, could be 
taken close to energies of N2p2 and B2pz states in re- 
spective atoms. However, for such choice, the value of 
Eg would be about 0.2 eV while in the boron nitride a 
wider energy gap (5.9 eV) is observed. 

This considerable separation between valence and con- 
duction bands is explained by difference in electroneg- 
ativities of boron and nitrogen which is responsible for 
a negative charge transfer in the direction of nitrogen. 
Then the covalent bond gains partially ionic character 
which produces a Madelung energy. By means of DFT 
calculations it has been found that for a freestanding h- 
BN sheet about 0.56 electrons are transferred from B to 
Ni^ This yields a sizable Madelung contribution to the 
lattice stability of 2.4 eV per 1x1 unit cell. The distin- 
guishability of boron and nitrogen leads to the insulator 
behaviour. The valence band is mainly constituted by 
the nitrogen sublattice while the conduction band by the 
boron sublattice. This property is refiected in a density 
of states obtained by the ah initio calculations/^ 

B. Bilayer graphene versus graphene on hexagonal 
boron nitride 

A secular problem for bilayer systems leads to a quartic 
equation, which can be solved by radicals. For given 
parameters one can provide the solutions explicitly as a 
fimction of A:. At K and K' points the condition 
holds, thus we can express the eigenvalues of for 
arbitrary parameters: 
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E4 = e^,l (22) 

In Fig. [2] we plot the band structures of the system 
described by the general hamiltonian (|lip for various 
choices of (0 and fixed hopping parameters. For sim- 
plicity we assume e"^ = and — ~e^^. We start 
with the freestanding BLG (Fig. [2](a), e"i = e''^ = 0). A 
difference between onsite energies (Fig. [2](b)) generates 
the 'Mexican hat' structure in the valence band which is 
also observed in bilayer graphene with substrate induced 
asymmetry."^- Further increasing of the asymmetry value 
leads to the touching of two valence bands at the K and 
K' points (Fig. [l](c)). For the values of larger than 
|7i|/'\/2 (Fig. [2](d)) the bands assigned to the eigen- 
values i?3 and Ei change their order, thus in the limit 
of graphene on h-BN a lower valence band gains N2pz 
character (Fig. [l](e)). 

The band structure of graphene/h-BN system pre- 
sented in Fig. [l](e) has been obtained for the values 
of e"^ = — chosen to assure a proper energy gap for 
the case of h-BN monolayer described in the previous 
section. In consequence, a gap between valence and con- 
duction bands is about 50 meV which is of the same order 
as the one obtained from DFT calculations for graphene 
on four layers of h-BN substrate.— In the vicinity of the 
Brillouin zone corners the band structure possesses nearly 
carbon character indicating a weak interaction between 
the layers. 

In the present considerations, we have shown that evo- 
lution from BLG towards the graphene on h-BN system 
is related to changes of the onsite asymmetry in the bot- 
tom layer. 

C. Bilayer system in an external electric field 

The electronic properties of a graphene bilayer in the 
presence of an external electric field applied perpendicu- 
larly to the system have been widely studied.®'^ We are 
interested here in the properties of bilayer system con- 
sisting of graphene and h-BN monolayers. The external 
perpendicular field gives rise to an electrostatic energy 
difference between these two layers which is parametrized 
by U. The inter layer asymmetry induced by an external 
field in the bilayer system shown in Fig. [T](c) can be 
described by the hamiltonian^ in the tight binding ap- 
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FIG. 2. (Color online) Band structure near K (or K') point for (a) graphene bilayer and (e) graphene on h-BN with e"^=- 
€^'^ = 2.95 eV. Graphs (b), (c) and (d) show the band structures for bilayer systems with different values of energy parameters 
e2p, and e2p, describing the evolution of band structure from graphene bilayer towards graphene on h-BN. The onsite energies 
sets used for these three graphs: (b) e"i = -e^i = 0.2 eV, (c) e"^ = -e'^i = |7i|/V2 and (d) e"i=-e^i = 0.32 eV. The hopping 
parameters of graphite 71 and 70 known from experimental datai^S fitted into the SWMcC model^'^ and 70 ~ 70 have been 
chosen for presented calculations. Blue dotted-dashed line denotes a band with the value E\ feQ ll9|) at the Dirac points; red 
dashed line is assigned to the value E2 (eg |20|) : olive dotted to the value £4 (ea l22|) and red solid to the value E3 (eg |2T|) . 




U[V] 



FIG. 3. (Color online) Variation of the energy gap as a function of energy difference between two layers. The insets show the 
band structures in three characteristic points indicated on the gap variation curve. The TB parameters are taken from Table U 
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FIG. 4. (Color online) Band structures of graphene on h-BN bilayer for different values of U. Blue dotted-dashed line denotes 
a band with the value Ei (ea l24|) at K point; red dashed line is assigned to the value E2 (eq. [25]); olive dotted to the value E4 
(eq. I27[l and red solid to the value E3 (eq[26]). The TB parameters are taken from Table ID 
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When the condition ([T5|) holds, one can receive the 
eigenvalues for arbitrary parameters: 
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where we have assumed, for simplicity, e 

As a starting point to study the influence of the elec- 
tric fleld on the band structure we flnd the values of pa- 
rameters U for which two bands touches in the K [K') 
point. A simple analysis shows that it occurs in three 



cases. First one gives the value of U when the bands 
corresponding to Ei and are joined: 



Ua 



(28) 



second one when the bands corresponding to E^ and -E4 
meet in K point: 
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and the third one comes from the condition E2—E4 
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Figure ini shows the changes of a band gap as a function 
of U while the insets present the band structure around 
the K point in three cases described by conditions (|28l) . 
(I29l) . ((50)) and indicated by letters A, B and C on the 
plot. Since the condition ([28]) gives an imaginary value 
of Ua , the band structure in point A on the curve can be 
considered only by taking its real part. Inset A in Fig. [3] 
depicts the band structure determined for Ua — — e"^ 
and it follows that it is the lowest value of the voltage, 
for which a direct band gap is observed. It achieves a 
value of 71 . The second condition ([29]) predicts that two 
valence bands touches, which is shown in inset B of Fig. [3] 
Finally, condition ([50]) leads to the situation when the 
valence and conduction bands meet in the corner of the 
Brillouin zone (inset C in Fig. [3]). Moreover, in the last 
case we observe a large slope of conduction band in the 
vicinity of K point which indicates that nearly conical 
touching can be recognized (m* = 3.4 • 
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In Fig. |4]we present the band structures for the values 
of voltages chosen in the way to show what happen in 
the regions outside the boundaries determined by — e"^ 
and — e^^ and in the regions between points A - B and 
B - C (see Fig. [3]). For the values of voltages lower 
than -e"i (Fig. |4](a)) and higher than -e^^ (Fig. |4](d)) 
the band gap is indirect which is not observed in case of 
the bilayer graphene when valence and conduction bands 
have the same Mexican hat shape. Moreover, for higher 
positive voltage (Fig. 111(d)) we observe the restructur- 
ing of conduction band which seems to be a result of an 
anticrossing interaction between highly localized states 
of electronegative nitrogen atoms with the extended con- 
duction band states of graphenei^ The electronic states 
loose their localized character and form band through an 
interaction with extended states. For the voltages from 
A - B range of values (Fig. |4](b)), we observe an up- 
ward movement of the valence bands edges. It can be 
interpreted as induced by hybridization of the carbon pz 
states comprising the higher valence band with the close 
lying localized Pz states of N. This bands anticrossing 
phenomenon leads to the shift of the valence bands edges 
and to the reduction of the band gap.^^ 

It is worthwhile to notice that the band structure mod- 
ifications of the bilayer graphene/h-BN system induced 
by the external electric field is similar to the behaviour 
of the band structure in mismatched III-N-V and III-V 
semiconductor alloys i^"— 
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III. DENSITY FUNCTIONAL CALCULATIONS 



We have performed density functional calculations for 
bilayer graphene/h-BN system using PWscf package^- of 
the Quantum Espresso distribution^. A plane wave ex- 
pansion with an energy cutoff at 400 Ry and the Perdew- 
Zunger LDA exchange - correlation functional"*?, were 
used. Core electrons were described by norm conserv- 
ing pseudopotentials. The Brillouin zone was sampled 
on a uniform 36 x 36 x 1 grid with the tetrahedron 
method^"^. The interactions with spurious replicas along 
the z direction were avoided by setting a large size of 
a supercell z ~ 16 A. We have taken a value of lattice 
constant a=2.456A. Total energies were convergent to 
10"-* Ry. 

We fit the TB band structure with DFT data to test 
the model analyzed in Sec. II B. The fitting procedure has 
been applied to the points taken from the straight vicinity 



TABLE I. Tight binding parameters calculated from best fit 
to DFT data. All values are given in electronvolts. The en- 
ergies ejp^ and e^^^ are measured from the zero level fixed by 
the condition e?ij = 0. 
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3.34 
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FIG. 5. (Color onfine) The TB fit along TKM line in the 
vicinity of K point for LDA calculations. Red dashed lines 
denote DFT bands while black dotted lines denote TB fit. 
The inset show DFT band structure in FK and KM direc- 
tions. The energies are measured from the zero level fixed by 
the condition e?^ = 0. 



of K point along TK line using the least squares method. 
We present the calculated DFT and TB band structures 
in Fig. [5l It is easy to see, that both electronic structures 
are well fitted. The small discrepancies start to appear 
only at the higher conduction band. The obtained values 
of TB parameters are listed in Table HI 

As it is expected, in local density approximation 
(LDA), the difference \Ei — £4] between boron nitride 
bands is underestimated. Its value should be higher than 
5.9 eV*"*, while the currently determined value is in the 
order of 4.7 eV which is in a good agreement with LDA 
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FIG. 6. DFT band structure along FKM line in momentum 
space of bilayer graphene/h-BN in the presence of external 
electric field with the effective value of U = 1.40 V applied in 
the direction perpendicular to the layers. The inset presents 
zoom of the band structure in the vicinity of K point, where 
the band gap is closed. 

results for pure boron nitride honeycomb^. 

In the second stage, we have performed the ab initio 
calculations for graphene/h-BN system in the presence 
of a constant electric field simulated by a sawlike po- 
tential. The results confirm TB predictions related to 
the band structure modifications described in Sec. lie. 
It should be noticed that similarly to the case of bi- 
layer graphene^^''''^ we observe screening effects. It means 
that in order to obtain expected shift and deformation 
of bands we have to apply an electric field higher than 
estimated in TB model. Figure [6] presents DFT band 
structure for the case given in inset C of Fig. |31 The TB 
predicted value of [/ = — e^^ is too low to reproduce the 
gap closing, in fact we have to apply U = 4.15 V. 

In Fig. [7] we show the band structure of the graphene 
on h-BN in the presence of an electric field parametrized 
by negative U of the effective value equal to -2 V. We can 
observe a direct gap of over 127 meV between valence and 
conduction band. One can notice, however, that at the 
r point one of the boron nitride originating bands starts 
to coincide with the energy level of the bottom edge of 
carbon conduction band (at the K point), which leads to 
an indirect band gap of the system. For the higher abso- 
lute values of potential difference, a band gap decreases 
due to the presence of unwanted edge and achieves nega- 




FIG. 7. (Color online) DFT band structure along FKM line 
in momentum space of bilayer graphene/h-BN in the pres- 
ence of external electric field with the effective value of U = 
-2.0 V applied in the direction perpendicular to the layers. 
Red arrow points out the edge of conduction band originat- 
ing from hybridized states comprising tt carbon band and a 
boron nitride band. 

five values typical for semimetals. When a system gains 
semimetallic character, very strong screening effects pre- 
vent further increase of the effective voltage. 

IV. FINAL REMARKS 

We have studied the electronic properties of a 
graphene/h-BN bilayer system using the minimal tight 
binding model. Particular attention has been given to 
the presence of external electric field perpendicular to 
the bilayer system which gives rise to a finite gap in the 
spectrum, whose opening and width are controlled by 
the applied voltage. The use of TB approach allows to 
show that graphene/h-BN system and bilayer graphene 
are very closely related. 

Moreover, we propose a very simple model as a start- 
ing point for similar considerations of more sophisticated 
substrates and interactions with metallic contacts4^i^ 
Further generalization on a larger number of AB stacked 
layers of graphene and AA stacked h-BN sheets can also 
be made. 

We have shown that the band structure may be mod- 
ified and tuned to meet specific requirements for poten- 
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tial logic devices. We report the energy gap tuning in 
the range starting from to about 130 meV and effec- 
tive masses in conduction band varying from 3.4- 10^'^ rric 
to 1.1 • 10~^mo, which are smaller than Dirac fermions' 
masses in bilayer graphene. One can easily predict that 
additional quasi-particle interactions included in calcula- 
tions, for instance within a GW scheme, will indicate an 
increased value of the gap.^^ Moreover, DFT as a ground- 
state theory does not give a proper description of excited 
states. Therefore, it is possible that an energy of bottom 
edge of the conduction band (see Fig. [7]) is not correctly 
determined. On the other hand, if one could engineer 
the band structure to shift or deform this conduction 
band edge then the voltage induced band gap could be 
as large as predicted by TB calculations (over 430 meV 
see Fig. 111(a)). 

Since a significant step towards the bilayer graphene 
FET production has been recently taken,'' it has to be 
stressed that various modified bilayer systems could play 
a crucial role in the future nanoelectronics. Recent stud- 
ies on transferring exfoliated graphene onto arbitrary 



substrates^ suggest that the fabrication of graphene 
on h-BN could be now realized. 
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